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O ' ABSTRACT 

OA \ In the very near future, weak lensing surveys will map the projected density of the 

universe in an unbiased way over large regions of the sky. In order to interpret the 
results of studies it is helpful to develop an understanding of the errors associated with 
quantities extracted from the observations. In a generalization of one of our earlier 
works, we present estimators of the cumulants and cumulant correlators of the weak 
lensing convergence field, and compute the variance associated with these estimators. 
, By restricting ourselves to so-called compensated filters we are able to derive quite 

J> ■ simple expressions for the errors on these estimates. We also separate contributions 

from cosmic variance, shot noise and intrinsic ellipticity of the source galaxies. 
00 . 
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■ 1 INTRODUCTION 

f The recent detections of weak gravitational lensing of background galaxy images by large-scale structure (Bacon et al. 2000; van 
JL , Waerbeke et al. 2000) provide an added impetus to the development of statistical methods for handling high-quality data from 
» l ■ weak lensing surveys which should provide us with valuable information about the mass distribution in the universe (Mellier 
' 1999; Bernardeau 1999; Bartelmann & Schneider 2001). The particular benefit of weak lensing surveys is that they permit 
, us to probe the distribution of underlying mass in a fashion that does not depend on an understanding of the relationship 
T* • between galaxies and the distribution of dark matter. 
, ' Following the directions set in earlier work by Gunn (1967), Blandford et al. (1991), Miralda-Escude (1991) and Kaiser 
, (1992), most current progress in weak lensing can broadly be divided into two categories. Villumsen (1996), Stebbins (1996), 
5_l ' Bernardeau et al. (1997) and Kaiser (1987) have focussed on the linear and quasi- linear regime by assuming a large smoothing 
angle. Several other authors have developed computational techniques to simulate weak lensing catalogs: numerical simulations 
of weak lensing typically employ N-body simulations, through which ray tracing experiments are conducted (Schneider & Weiss 
1988; Jarosszn'ski et al. 1990; Lee & Paczyn'ski 1990; Jarosszn'ski 1991; Babul & Lee 1991; Bartelmann & Schneider 1991: 
Blandford et al. 1991). Building on the earlier work of Wambsganns et al. (1995, 1997, 1998) the most detailed numerical 
studies of lensing have been made by Wambsganns, Cen & Ostriker (1998). Other recent studies using ray-tracing experiments 
have been conducted by Premadi, Martel & Matzner (1998), van Waerbeke, Bernardeau & Mellier (1998), Bartelmann et al. 
(1998), Couchman, Barber & Thomas (1998) and Jain, Seljak & White (2000). While the former, perturbative analyses can 
provide valuable information at large smoothing angle, this approach can not be used to study lensing on small angular scales, 
because the perturbation series involved start to diverge. 

More recent studies (Hui 1999; Munshi & Jain 2000, 2001; Munshi & Coles 1999, 2000; Valageas 1999a,b) have demon- 
strated that, in the highly non-linear regime, it is possible to combine the well-motivated hierarchical ansatz (Davis & Peebles 
1977; Peebles 1980; Fry 1984; Fry & Peebles 1978; Szapudi & Szalay 1993, 1997; Colombi et al. 1997; Scoccimarro & Frieman 
1998; Scoccimarro et al. 1998; Balian & Schaeffer 1989; Bernardeau & Schaeffer 1992; Bernardeau & Schaeffer 1999) with 
the scaling relation for evolution of two-point correlation functions (Hamilton et al. 1991; Nityananda & Padmanabhan 1994; 
Jain, Mo & White 1995; Padmanabhan et al. 1996; Peacock & Dodds 1996) to make very accurate predictions of the statistics 
of the lensing convergence field for very small smoothing angular scales. This approach offers an analytic route to the study 
of small scale lensing to complement the numerical approaches mentioned above. 
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The hierarchical ansatz has been used to show that lower-order moments such as cumulants and cumulant correlators can 
be modelled very accurately. It has also been found that the probability distribution function (PDF) and the bias associated 
with hot spots in convergence maps can also be predicted very accurately using this formalism. Higher order moments are 
more sensitive to the tail of the distribution function they represent, and are consequently more sensitive to measurement 
errors arising from finite size of the sample. Although there have been many detailed studies to quantify measurement errors 
for moments of the density field (Colombi et al. 1995; Colombi et al. 1996; Szapudi & Colombi 1996; Hui & Gaztanaga 1998) 
similar studies for weak lensing surveys are still lacking. Schneider et al. (1998) proposed different estimators for extracting 
the variance from convergence maps and the errors associated with them. In this paper we extend such results to incorporate 
all higher order cumulants and their two-point counterparts the cumulant correlators (Szapudi & Szalay 1997) which were 
used in the context of lensing by Munshi & Coles (2000). We study the contribution to the error involved in using these 
estimators by computing their variance. We list contributions from different sources (including the discrete nature of the 
source distribution, the intrinsic ellipticities associated with source galaxies, and the finite size of the catalogue). In previous 
studies of error estimations the higher order cumulants were assumed to be zero, as there has not been until recently an 
analytic prediction for the hierarchical parameters Sn in the highly non-linear regime. Combining our result with recent 
analytical prediction for Sn parameters for small smoothing angles will provide an accurate way to compute estimation errors 
and hence actual possibility of measuring these quantities from observational data. 

The layout of the paper is follows. In Section 2 we introduce the estimators, and in Section 3 we explain the different 
types of averaging involved in computing the dispersion and mean of these estimators. In Section 4 we develop a diagrammatic 
formalism to compute the mean and the dispersion and derive very general expression for scatter in estimates of the cumulant 
correlators of arbitrary order for an arbitrary number of points. In Section 5 we discuss the importance of our results in a 
general cosmological context. We have presented the detailed expressions for specific lower-order moments in an appendix for 
easy reference. 



2 ESTIMATORS FOR CUMULANTS AND CUMULANT CORRELATORS 

The statistics most frequently used to quantify the nature of clustering from galaxy catalogs are the moments of various 
orders. These are useful both to quantify the nature of non-Gaussianity and also to constrain the nature of initial conditions. 
A particularly useful way of combining moments is in the form of cumulants, which have been used to quantify both galaxy 
clustering and lensing surveys. Unlike the cumulants derived from galaxy catalogs, cumulants of lensing fluctuations can also 
differentiate between different cosmological models. These will be the most useful statistical descriptors for future weak lensing 
surveys. 

Schneider et al. (1998) have proposed the use of aperture mass statistics based on the use of a compensated filter function 
U to smooth the weak lensing convergence field k defined over a circular patch of sky with a radius 6q\ 

M ap (6> ) = J d 2 9U{9)n{9). (1) 

This particular filter function has many useful properties which allows us write M ap in terms of the measured tangential 
component of the shear jt inside a circle of radius 9 on the sky: 

M ap (0 o ) = / d 2 9Q{8) lt {9), (2) 



where Q{9q) and U(8) are related by 

= T2 f d9'U(9') - U(9 ) (3) 
8o J 

(Schneider et al. 1998). We will use the second definition of M ap in our analysis. The analysis can be extended to any other 
specific form of window function, such as a top hat, although it will no longer possible to directly relate the smoothed 
convergence fields with the galaxy shear 74 measured observationally. 

We begin by defining an estimator for the cumulants which is a natural generalization of the lower-order estimators used 
by Schneider et al. (1998): 
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where n is the number of galaxies in the patch of size tv9q, the e; are individual image ellipticities, and N is the order of the 
cumulant. The function Q and its relation to the compensated filter function U are defined above. In equation (4) we use the 
notation 
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We propose a family of new estimators for cumulant correlators. Cumulant correlators were introduced in the context of 
galaxy surveys by Szapudi & Szalay (1997). The new estimators are defined as: 
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This approach can in principle be extended to s-point cumulant correlators which are defined over s different patches of 
the sky where measurements have been conducted. 
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For detailed description of these quantities see Munshi et al. (2000) context of galaxy surveys and Munshi & Coles (2000) for 
their weak lensing counterparts. It is well-known that these quantities carry more information then their one point counterparts 
the cumulants. 

In order to be useful, the signal-to-noise ratio involved in measurements of these quantities should be high. Our main aim 
in this paper is to develop analytical results which take into account contributions from various sources of error (or "noise"). 
These include the distribution of intrinsic ellipticities for the lensed galaxies, shot-noise resulting from the discrete nature of 
galaxy distributions, and the finite size of the catalogues. While the last two contributions also arise during the analysis of 
projected galaxy catalogs, the intrinsic ellipticity distribution is a source of uncertainty unique to weak lensing surveys. 

Finite catalogue size has two principal effects. The first is that a finite volume (or projected area) can not reveal information 
about structures on a scale larger than the sample. In cosmic microwave background studies this is sometimes, though not 
entirely accurately, referred to as "cosmic variance". We can generally neglect this first aspect of finite volume by using 
compensated filters. The other main effect is more difficult and involves the effect of sample boundaries or edges imposed by 
sample geometry, particularly if it involves complicated masks. We shall generally refer to this second source of error as "finite 
volume" . 

We will show that the different contributions to measurement errors of cumulant correlators are in general factorizable 
and can be separated into "pure" terms and hybrid errors associated with measurements of one-point cumulants. Our results 
are quite general and are valid for arbitrary order and for arbitrary number of smoothed patches. Using the rules we have 
developed in this paper it will also be possible to compute the higher order moments of errors associated with different 
estimators. For related discussions, see Szapudi & Colombi (1996) and Szapudi, Colombi & Bernardeau (1999). 



3 DIFFERENT TYPES OF AVERAGING 

The total shear et can decomposed into the part which is due to the intrinsic source ellipticity e^ and the ellipticity introduced 
by distortion of images due to weak lensing, which we will denote by 7t. To proceed further we need to consider the effects of 
three different types of averaging process. First, there is the average over positions of source galaxies within the patches 
which are used to compute the cumulant correlators. Let us denote this operator, operating on an arbitrary statistic M, by 
P(M) where 

w=(n/0)«- w 

Second, we have the average over the distribution of intrinsic ellipticities. This particular source of noise is generally 
assumed to be Gaussian. It is further assumed that the ellipticities of neighbouring galaxies (in projection) do not correlate 
with each other. For a given estimator say M we will denote this average by G(M). It will only operate on the intrinsic 
ellipticity variables i.e. e s ti . So we can write: 

a 1 a 2 
G(euetj) = Sij-^-, G(e u etjetk) = 0, G(€ue t ]e tk € tl ) = 5ij5 k i -y + cycl.perm.; G(7i7j . . .) = 7»7j ■ ■ ■ ( 9 ) 

Finally, there is the ensemble average over different realizations of the sky. The ensemble averaging is commonly denoted 
by (...}. Putting these averages together we can write the expectation value of a given statistic M for a particular patch on 
the sky as 

E(M) = (G(P(M))). (10) 



For a more detailed description of these operators and their commuting properties see Schneider et al. (1998). 



Figure 1. Computation of the variance by a diagrammatic technique. Each circle represents one copy of the smoothed patch of the sky 
in a particular direction. To compute the variance we need to consider two copies of each patch, and so on to higher order. The diagram 
shown above represents one particular term in the expression for the dispersion. The points in each patch represents one et, which can 
represent an intrinsic source ellipticity e a or contribution due to weak lensing 7. By assumption, source ellipticities in different patches 
do not correlate but the ellipticities induced by weak lensing do correlate. Each line joining 7 pairs has a weight factor of M s and the 
lines joining e's represent Mg. Points which are not joined by different lines denote M ap . Permutations of these terms in one patch of 
the sky do not depend on the other patches. We have shown an error term (without subtracting the disconnected parts) for two-point 
cumulant correlators C73 but it can be generalized trivially to incorporate multi-point smoothed cumulant correlators. 

4 MEAN AND DISPERSION OF ESTIMATORS 

The basic formalism we will adopt in our analysis of the bias and dispersion of these estimators is similar to that developed 
by Schneider et al. (1998). Throughout the following analysis we will focus exclusively upon compensated filters. It was shown 
by previous studies that such filters are local in nature, meaning that correlations between neighboring cells will be quite 
negligible. As we shall see, our analytical results for the appropriate error terms for the one-point cumulants contain only 
higher-order one-point cumulants, and not two-point cumulants. Although our results will be quite accurate for the specific 
case of compensated filter functions, for other smoothing windows it is necessary to consider the correlation between different 
smoothing windows. This introduces a bias in our estimators in addition to the scatter we are studying here: for a thorough 
discussion of the ramifications of this, see Bernardeau et al. (1997). 

We will consider several patches of the sky where measurements of shear are performed, and we will do an averaging 
over galaxy positions within each of these patches, intrinsic galaxy ellipticity and finally an ensemble averaging over all sky- 
positions. These results then will be generalized to the case when simultaneous measurements are carried over many s-tuplcs 
of patches for cumulant correlators of order s. 

To compute the mean or dispersion of these estimators we use a diagrammatic technique, which will simplify the compu- 
tation and will allow us to write a very general expression for the dispersion associated with them. 

4.1 Rules 

To compute the dispersion we have to consider an identical copy of the same patch. After expanding the multinomial expression 
that results from this splitting into intrinsic and lensing- induced shear we can express the statistic as a sum of various terms 
which are just products of various combination of powers of ef^ and 7 t represented as a diagram as shown in Figure 1. To 
compute the third order moment of errors it is necessary to consider three copies of the same patch, and so on. The action of 
the various operators as discussed above will result in pairing of these stochastic variables which can be computed by following 
the rules listed below. 

The different pairings of e will be considered between copies of the same patch. The different types of pairing and the 
rules for dealing with them are: 

• 7 pairing: each of these terms will contribute one M 2 term which will denote the contribution to the error the from 
discrete nature of the galaxy distribution (denoted by solid lines in Figure 1). Here M s is defined by 

M 2 = ^? 2 [ d 2 6 Q 2 {6) 1 2 t {6) (11) 

• e pairing: each of these pairings will contribute a M 2 term, which arises from the intrinsic ellipticity of galaxies (denoted 
by dashed lines in Figure 1). 

Mg = 7r#o f d 2 9 Q 2 (8). (12) 

The intrinsic ellipticity distribution is assumed to be Gaussian and uncorrelated with other patches, as mentioned above. 

• No pairing of 7 or e which will take contribution from correlation term between different patches (denoted by black 
dots in Figure 1). The amplitude associated with these terms are simply M ap which we have defined before. 

Total number of pairs which can be made out of n objects is n!: out of these pairs some of them will involve 7 pairing 
and some of them will correspond to e pairing. 




Figure 2. M ap for various cosmological models as a function of smoothing angle 9o which we have used in our error estimation. 
Parameters characterizing these cosmological models can be found in Munshi & Jain 2000. 

4.2 Expectation Values 

It is not difficult to show that the estimators we proposed in Section 2 for the cumulants and cumulant correlators are unbiased 
estimators. From a diagrammatic point of view, it is clear that there are no terms that correspond to coupling of e^ s ' terms 
with either from the same patch (because no indices in the same patch can be equal) or from two different patches (because 
they are not correlated). We also note that 7t and are also uncorrelated. This means that the only contribution comes 
from terms in which the y t are correlated. Hence we can write, for the cumulants, 



E(P(A(M N ))) = (f[ J^eiQ(0ih t {Oi)\ = (Mg) 



and similarly for cumulant correlators, 
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which shows that they are unbiased estimators of multi-point moments. 



4.3 Dispersion 

As we discussed above, for the computation of variances we have to consider two copies of the same patch. Following the rules 
given above we can finally write down the expression for the dispersion of an iV-th order cumulant as 
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The symbol C p is used to denote the number of combinations of p objects taken from N. The first term here represents 
the case when all points are 7 t in both patches and there is no pairing of these points within the copies of the same patch. 
The second term represents only 74 pairing between pairs of 74 from two different copies of the same patch. The third term 
represents the e t pairing within the copies of the same patch and the last terms is a mixture when some of the couplings are 
et coupling and some of them are 7 t coupling. 

A similar analysis can be performed for s-point smoothed cumulant correlators of order Ni + . . . + N s . We find 
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Figure 3. The parameter p(0o) which appears in our expressions for fractional errors of various cumulants. Different curved denote various 
background cosmology. Notations arc same as in previous figure. For small values of p(8o) the errors are dominated by contributions from 
shot-noisc where as for large values of p(0q) it is dominated by cosmic variance. Fractional errors in various parameter are polynomial in 
inverse p and minima of these polynomials occur on an arcminute scale. For left panel we have assumed 30 galaxies per square arcminute 
and 60 galaxies per square arcminute for the right panel. 
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The above expression is the general expression for the dispersion which is derive using the approximations we have explained. 
At each order, for a given family of cumulant correlators, the contribution to dispersion originates due to the three different 
effects explained before and their cross terms. 

Specific expressions for the dispersion of estimators for particular cumulants and cumulant correlators are given in the 
Appendix. Our results here are derived for the case when there is only one field of view. In practice, however, the part of the 
sky observed will be covered by several (possibly overlapping) fields of view from which the cumulant and cumulant correlators 
are estimated. For the case of compensated filter the generalization is quite straightforward but for other filters it will require 
a more rigorous analysis than is presented here. In the context of galaxy clustering, such effects have already been explored by 
Szapudi & Colombi (1996). The effect of source clustering will also have to be taken into account for more realistic estimation 
of errors. 



4.4 Errors associated with Sn parameters 

For practical purposes the normalised cumulants and cumulant correlators generally involve 

Sn <m! p >"-i • (17) 

These quantities are of great importance as they can quantify the non-Gaussianity in the weak lensing field and also can be 
used to study the effects of background cosmology on weak lensing statistics. Although the Sn parameters for weak lensing 
surveys have already been measured using top-hat filters, there has been not work in this direction for computation of Sn 
parameters using compensated filters for high N. 

Our error expression can be written in a compact form if we use the parameters Eat parameters which can be related to 
the Sn parameters by following expressions. These definition of Sjv also include the disconnected contributions. 

<M a 4 p > = 5 4 <A4 2 p) 3 + 3(M a 2 p ) 2 = S 4 <M a 2 p } 4/2 
<M a 5 p ) = S 5 {M!X + 10S 3 <M a 2 p } 3 = S 5 (M a 2 p ) 5/2 

<M a 6 p ) = S 6 (M! P ) 5 + 15S 4 (M a 2 p > 4 + 105 3 2 <M ap ) 6 = E 6 (M a 2 p ) 6/2 (18) 

Using the expressions derived in previous sections and in the limit when number of source galaxies are much higher compared 
to unity we can write down the following expressions for variance in M3 and M2. 
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Figure 4. Estimation error for computation of (M 2 p ) as a function of smoothing angle 9q. A survey area of 1° is assumed for left panel 
and a survey area of 8° is assumed in the right panel. Different curves denote different background cosmologies. Line style is the same as 
that of Figure 2. While measurement errors are dominated by the shot noise (which are induced due to intrinsic ellipticity distribution 
of the galaxy) and the Poisson shot noise for smaller smoothing angles. In large smoothing angles it is more dominated by finite size of 
the survey area. Estimation errors for the variance depends on the higher order moments. In particular we found that the normalized 
cumulant 54 appears. Unfortunately there are no theoretical estimates for computation of S4 using compensated filters. Earlier results 
in the literature used S4 = 0, which we have used here. Clearly the fractional errors are quite small even for moderate survey size and 
for large sky coverage power spectrum estimation can be carried out with very high precision. All our calculations are done under the 
hypothesis that the correlation length scale associated with compensated filter is quite small compared to length-scales associated with 
field of view. This was demonstrated by Schneider et al. (1998). For top-hat filters correlations between neighboring cells will have to be 
taken into account. In our computation the intrinsic eliipticity of the source galaxy cr e is set to be equal to .2. 




Figure 5. Estimation error for computation of S3. We have assumed that all sources are at a fixed red-shift. A survey area of radius 
1° (left panel) and 8° (right panel) are assumed. Number density of source galaxies n survC y in the left panel is 30 per arminute 2 in the 
right panel it is fixed at 60 per arcminute 2 . Error is plotted as a function of the smoothing angle 8q. Different curves denote different 
background cosmolgies. Line style is same as in previous figures. While the errors are mostly dominated by shot noise at a very small 
resolution angle, finite size of the survey are more pronounced for larger smoothing angles. However for smoothing angular scales of 
few arc minutes the fractional errors are quite small for reasonable sky coverage. This is more interesting as non-Gaussian estimators 
such as S3 can be used to distinguish between different cosmologies. Notice that fractional errors do not change much while we change 
the background cosmolgy. For computation of the fractional errors in S3 we have used the analytical results for compensated filters by 
Bcrnardcau & Valageas (2000) which are in agreement with numerical computations by Rcblinsky ct al.(1999). Although these results 
are for smoothing radius 80 = 4', the Sjv parameters are known not to change much at very small angular scales where we can assume 
them approximately constant for computations of error bars. 
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Where we have introduced the quantity p defined by the following expression: 
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For small p the errors are dominated by shot noise i.e. in the small smoothing angle regime, while for large p the errors are 
dominated by cosmic variance. The error terms are a polynomial function in 1/p. The minima of these polynomial occurs 
on a scale of a few arcminutes, where neither cosmic variance nor the shot noise terms are dominant for reasonable survey 
strategy. 

The behaviour of M ap for various cosmologies is shown in Figure 2. The function p is plotted in Figure 3. Using the 
results we have obtained, the scatter associated with the Sat parameters can now be written in terms of these quantities: 

Sn - <M£> + {N 1] (M2 P ) 

For multiple fields covering the whole field of view factors of y/~N f (where Nf is the number of cells covering the field) will 
have to be incorporated (see Schneider et al. 1998). For example, the variance is a factor Nf smaller. This, of course, is only 
the case for compensated filters like those we are studying here. Effect of over-sampling has not been studied in the context 
of weak lensing and needs a dedicated analysis. 

Final results for the scatter in variance and S3 are plotted in Figure 4 and Figure 5 respectively. For specific case studies 
we have assumed two different hypothetical survey strategies. In one case we assume that the radius of the survey is 1° in 
the second case it is fixed at 8°. The number density of source galaxies in the first case is assumed to be 30 arcmin~ 2 and 
twice that value in the second case. The parameter a e which characterizes the intrinsic ellipiticity distribution of galaxies is 
fixed at .2 in both cases. While the first case is representative of various ground based surveys with less exposure time which 
are currently in progress, the second case is more idealistic and hope to be achieved by future weak lensing survey with long 
exposure time probably from outside earths atmosphere. In both cases the errors are minimum for arcminute scale. Our study 
shows that skewness parameters S3 can be measured very accurately from second type of surveys although signals might be 
detectable even in first case with considerable uncertainty resulting in large error bars. This underlines the importance of 
future weak lensing surveys from cosmological perspective. 

Direct determination of higher order Sat parameters for the case of weak lensing surveys has so far only been reported 
by Bernardeau et al. (2002). Higher moments become more difficult to measure because better statistics are needed. As the 
expression for the error in such quantities contain Sn parameters of higher order we need to make certain order-of-magnitude 
approximations to make progess in estimating errors. Schneider et al. (1998) assumed that E4 = 3. In addition to that we 
have assumed that Eg = E§ for error expression in S3 . A more accurate error estimation will have to wait till we can compute 
these quantities analytically. However any analytical results will have to be checked with numerical result for which we need 
to have a better understanding of their error properties. Our results are only a first step in this cycle and will have to be 
refined by future studies. 

We stress again that our derivation is only exact for the case of compensated filters when neighboring smoothed regions 
are not correlated. For the case of other filter functions such as the top-hat filter a more general treatment is necessary. We 
plan to present such results elsewhere. We have also assumed that the source redshifts are determined with high accuracy. 
Lack of knowledge of source redshifts is a potential source of additional error which too needs to be taken into account for 
more realistic estimation of cosmological parameters from future weak lensing surveys. 



5 DISCUSSION 

It has been the purpose of this paper to estimate the errors involved with the extraction of statistical information from 
weak lensing surveys. We have focussed on the cumulants and cumulants correlators, which are normalized moments of the 
smoothed one- and two-point probability distribution functions. These quantities are widely used to quantify the statistical 
nature of clustering of the mass distribution in the study of galaxy surveys. In this context, estimators of these statistics are 
prone to error from finite catalogue size and Poisson (discreteness) effects. The application of similar methods to weak lensing 
studies is clearly appropriate, but introduces an additional source of error. This paper allows for these additional error terms. 

Various estimators have already been proposed for the computation of the variance and skewness of cosmic shear smoothed 
with a particular window function. We have also generalized these suggestions to statistics of arbitrary order via a generalized 
estimator which estimates the statistics of smoothed convergence field in an unbiased way. We have also proposed a new set 
of estimators which are useful for measuring cumulant correlators, and which are natural generalizations of their one-point 
counterparts. We have shown that our estimators constitute a family of unbiased estimators for cumulant correlators. 

We have also computed the dispersion of these estimators, which is essential to determine the signal-to-noise ratio 
associated with them. We have found compact expressions for the dispersions for arbitrary order and also for arbitrary 
number of points. We have been also able to separate contributions from different sources of noises such as the finite size of 
the galaxy catalog, finite width of ellipticity distribution of source galaxies and Poisson noise due to finite number of source 
galaxies in the field of view. Our results are valid for both large and small smoothing angles. In case of large smoothing 



angles, where perturbative calculations are still valid, one need to use the quasi-linear values of Sn and C pq parameters 
associated with the expressions for finite volume corrections. For smaller smoothing angles we have to replace these numbers 
by forms suitable for the highly non-linear regime, which have already been computed by several authors recently based on 
the hierarchical ansatz (Hui 1999; Munshi & Jain 2000, 2001; Munshi 2000, Valageas 2000a,b). 

The effect of source clustering, which we have ignored, will also introduce corrections terms in the measurement of lower 
order moments. We have also ignored the effect due to lens coupling. Some of these issues have been studied by Bernardeau 
et al. (1997) and Schneider et al. (1998) with the result that such corrections seem to be negligible, at least in the quasi- 
linear regime. Studying the effects of source clustering using ray-tracing experiments in the highly non-linear regime is difficult 
because most such simulations propagate light rays backward: the source position is consequently left arbitrary and determined 
only by lensing due to the intervening mass. 

The validity of the Born approximation, which underpins lensing calculations, has also been studied in the quasi-linear 
regime. It has been shown that corrections arising from higher-order terms in the photon propagation equation are negligible 
in quasi-linear regime. Similar conclusions have also been found to be valid in the highly non-linear regime by comparing 
ray-tracing simulation against analytical results using hierarchical ansatz. (Hui 1999; Munshi & Jain 2000, 2001; Munshi 2000; 
Valageas 2000a,b). 

We have also assumed that the galaxy intrinsic ellipticities are not correlated but it may be possible that the galaxies 
are not randomly oriented and there may be a a coherent alignment due to the geometry of the large-scale structure in which 
they are embedded (e.g. Heavens, Refregier & Heymans 2000). So far however no convincing observations of nearby structures 
have indicated that such an alignment exists (e.g. Mellier 1999) although several attempts have been made to unearth one. 



ACKNOWLEDGMENT 

DM was supported by a Humboldt Fellowship at the Max Planck Institut fur Astrophysik when this work was performed. It is 
a pleasure for DM to thank Bhuvnesh Jain, Francis Bernardeau, Patric Valageas. Katrin Reblinsky and Alexandre Refregier 
for many helpful discussions. 



REFERENCES 

Babul A., Lee M.H., 1991, MNRAS, 250, 407 

Bacon D.J., Refregier A., Ellis R.S., 2000, MNRAS, 318, 625 

Balian R., Schaeffer R., 1989, A& A, 220, 1 

Bartelmann M., Huss H., Colberg J.M., Jenkins A., Pearcc F.R., 1998, A&A, 330, 1 
Bartelmann M., Schneider P., 2001, Phys. Rep., 340, 292 

Bernardeau F., 1999, in Lachiczc-Rey, ed., Proc. Cargese Summer School on Theoretical and Observational Cosmology. Kluwcr, Dordrecht. 

Bernardeau F., Schaeffer R., 1992, A& A, 255, 1 

Bernardeau F., Schaeffer R., 1999, A& A, 349, 697 

Bernardeau F., van Waerbeke L., Mellier Y., 1997, A& A, 322, 1 

Bernardeau F., van Waerbeke L., Mellier Y., 2002, A& A, 389, 28 

Blandford R.D., Saust A.B., Brainerd T.G., Villumscn J.V., 1991, MNRAS, 251, 600 

Boschan P., Szapudi I., Szalay A.S., 1994, ApJS, 93, 65 

Colombi S., Bernardeau F., Bouchet F.R., Hernquist L., 1997, MNRAS, 287, 241 

Colombi S., Bouchet F.R., Schaeffer, R., 1995 ApJS, 96, 401 

Colombi S., Bouchet F.R., Hernquist L., 1996, ApJ, 465, 14 

Couchman H.M.P., Barber A.J., Thomas P.A., 1998, MNRAS, 308, 180 

Davis M., Peebles P.J.E., 1977, ApJS, 34, 425 

Fry, J.N., 1984, ApJ, 279, 499 

Fry, J.N., & Peebles, P.J.E., 1978, ApJ, 221, 19 

Groth, E., & Peebles, P.J.E., 1977, ApJ, 217, 385 

Gunn, J.E., 1967, ApJ, 147, 61 

Hamilton A.J.S., Kumar, P., Lu, E. & Matthews, A., 1991, ApJ, 374, LI 

Heavens A.F., Refregier A., Heymans C, 2000, MNRAS, 319, 649 

Hui L., 1999, ApJ, 519, 622 

Hui L., Gaztanaga E., 1999, ApJ, 519, 622 

Jain B., Mo H.J., White S.D.M., 1995, MNRAS, 276, L25 

Jain B., Scljak U., 1997, ApJ, 484, 560 

Jain B., Van Waerbeke L., 2000, ApJ, 530, LI 

Jain B., Seljak U., White S.D.M., 2000, ApJ, 530, 547 

Jaroszyn'ski M., Park C, Paczynski B., Gott J.R., 1990, ApJ, 365, 22 

Jaroszyn'ski M. 1991, MNRAS, 249, 430 

Kaiser N., 1992, ApJ, 388, 272 

Kaiser N., 1998, ApJ, 498, 26 

Limber D.N., 1954, ApJ, 119, 665 

Lee M.H., Paczyn'ski B., 1990, ApJ, 357, 32 

Miralda-Escude J., 1991, ApJ, 380, 1 



Melott A.L, 1998, preprint, astro-ph/9801011 
Coles P., Melott A.L., 1999a. 
Coles P., Melott A.L., 1999b 
Melott A.L., Coles P., 1999c. 



307, 
310, 
311, 



387 
892 
149 



Mo H., White S.D.M., 1996, MNRAS, 282, 347 

Munshi D., Bernardeau F., Melott A.L. , Schaeffer R., 199 9. MNRAS, 30:;. m 
Munshi D. 
Munshi D. 
Munshi D. 
Munshi D. 
Munshi D. 
Munshi D. 
Munshi D. 
Munshi D. 
Munshi D. 
Mellier Y.. 



Coles P., 
Coles P., 
Jain B. 
Jain B. 
Jain B. 



MNRAS, 
MNRAS, 
MNRAS, 
, 148 
, 797 



2000, MNRAS, 313 
2002, MNRAS, 329 

2000, MNRAS, 318, 109 

2001, MNRAS, 322, 107 
2000, MNRAS, in preparation 

1999, ARA& A, 37, 127 
Nityananda R., Padamanabhan T., 1994, MNRAS. 
Padmanabhan T., Cen R., Ostriker J. P., Summers. 
Peebles P.J.E 
Peackock J. A 



271, 976 

F.J., 1996, ApJ, 466, 604 
1980, The Large Scale Structure of the Universe. Princeton University Press, Princeton 
Dodds S.J., 1996, MNRAS, 280, L19 



Prcmadi P., Martel H., Matzner R., 1998, ApJ, 493, 10 
Press W.H., Schechter P. 1974, ApJ, 187, 425 

Reblinsky K., Kruse C, Jain B., Schneider P., 1999, A& A, 351, 815 
Schneider P., Weiss A., 1988, ApJ, 330,1 

Schneider P., van Waerbeke L., Jain B., Kruse G., 1998, MNRAS, 296, 873, 873 

Scoccimarro R., Colombi S., Fry J.N., Frieman J.A., Hivon E., Melott A.L., 1998, ApJ, 496, 586 

Scoccimarro R., Frieman J., 1998. ApJ. 520. 3 5 

Stebbins A., 1996, preprint, istro-ph/9609149 

Szapudi I., Colombi S., 1996, ApJ, 470, 131 

Szapudi I., Colombi S., Bernardeau F., 1999, MNRAS, 310, 428 

Szapudi I., Colombi S., Jenkins A., Colberg J., 2000, MNRAS, 313, 725 

Szapudi I., Szalay A.S., 1993, ApJ, 408, 43 

Szapudi I., Szalay A.S., 1997, ApJ, 481, LI 

Szapudi I., Szalay A.S., Boschan P., 1992, ApJ, 390, 350 

Valageas P., 2000a, A& A, 354, 767 

Valageas P., 2000b, A& A, 356, 771 

van Waerbeke L., Bernardeau F., Mellier Y., 1999, A& A, 342, 15 

van Waerbeke, L., et al. 2000, A& A, 358, 30 

Villumsen J.V., 1996, MNRAS, 281, 369 

Wambsganss J., Cen R., Ostriker J. P., 1998, ApJ, 494, 298 

Wambsganss J., Cen R., Xu C, Ostriker J.P., 1997, ApJ, 494, 29 

Wambsganss J., Cen R., Ostriker J. P., Turner E.L., 1995, Science, 268, 274 

Wang Y., 1999, ApJ, 525, 651 



APPENDIX A: EXPLICIT EXPRESSIONS FOR ERROR TERMS 



Al Cumulants 

Error terms for the one-point smoothed cumulants involve three contributions. For example, for the lowest order cumulant, 



cr (Ml) 



n(n - l)(M a 2 p ) + n{M 2 s ) + nM 2 g -± 



<M ap } 2 



(Al) 



The first of these contributions derives from finite volume effects included in the higher order moments such as M ap . For 
example the dispersion in variance depends on the fourth moment, and the dispersion in skewness will similarly depend on 
the sixth order moment. The terms denoted by M 2 are related to the fact that we have a finite number of galaxies: this will 
vanish if we take the limiting case of infinite number of galaxies. The terms with Mg are due to the finite width of the intrinsic 
galactic ellipticity distribution and would vanish in the limit in which there are an infinite number of galaxies in the survey. 

For the higher-order cumulants there are differences. Unlike the expression for the dispersion of Mi, we will have mixed 
terms in the expression for the dispersion in M2 which are denoted by various products of M s , M ap and M g . The following 
expression was derived by Schneider et al. (1998): 
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At third order we get the following expression for dispersion of A/3 
1 
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Notice that we can write (Fry 1984) (Mf p ) = S 6 (M% p ) 5 c + 155 4 (Mf p )* + lOSf (M 2 p ) 4 + 15(M 2 p ) 3 . For top-hat window functions 
we have analytic expressions for all lower order Sn parameters which can be used to estimate the effects of finite volume in 
the highly non-linear regime. For other window functions there is no such analytical expressions. 

Generally speaking, for large values of n, only the dominant contributions are considered and the rest are neglected. On 
the smaller angular scales the finite width of the galaxy ellipticity distribution dominates, and for very large smoothing angles 
it is the finite volume effect that dominates. Hence one can often neglect the terms containing M s altogether. 

At fourth order we get 

a 2 {Mi) = 1 
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In the statistical study of large scale distribution of galaxies generally cumulants are normalized by dividing them with 
suitable power of two-point cumulant or by the variance, for example in the construction of the Sn parameters. Since both the 
numerator and denominator are both affected by errors we have discussed above they will introduce a ratio bias as discussed 
by Hui & Gaztanaga (1999). 

Finally in this section we mention that the above results do not depend on the scale of non-linearity probed by weak 
lensing, but the appropriate values of Sn must be used. For example, when large smoothing angles are considered we should 
use the quasi-linear values of Sn parameters of the convergence map. 



A2 Cumulant Correlators 

To compute the cumulant correlators we consider two patches in the sky in the direction of 71 and 72. The results are very 
similar to the case of one-point smoothed cumulants. As before the errors associated with measurements can be differentiated 
in three types, i.e. the finite volume effects, errors due to intrinsic ellipticity of the source galaxies and errors associated with 
finite number of galaxies. The lowest order cumulant correlator is of course the smoothed two-point correlation function. The 
error associated with its measurement can be expressed as 

ct 2 (Mii) = — [{ni(ru - l)n 2 (n 2 - l)<M a 2 p ( 7 i)M 2 p (72)) + nin 2 (n 2 - l)(M 2 ( 7 i)M a 2 p ( 7 2)) 
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We have taken account that the number of galaxies in two different patches can be different. Notice that error now contains 
terms which are mainly the correlation of measurement errors in two different patches. 



Terms of higher order begin with M21 , which can be given in the form 
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These expressions are cumbersome, but worth listing because these are the cumulant correlators most likely to be mea- 
surable. In the above expressions we used 



In deriving the above expressions for measurement errors in two-point cumulant correlators we have again assumed that 
the intrinsic ellipticities of galaxies do not cross correlate among different patches. At increasingly higher orders the expressions 
for error contribution becomes more complicated, although as explained above for large number of galaxies we can take the 
m — ► 00 and ri2 — > 00 limit which will simplify these expressions. 

As in the case of cumulants normalized cumulant correlators can be derived by dividing cumulant correlators by suitable 
powers of variance within these patches and the correlation among these patches (see Munshi & Coles 2000). Normalized 
cumulant correlators are also denoted by C vq and in addition to the errors we have already discussed we will have ratio bias 
too. 

The formalism we have developed above can also be extended to compute the skewness associated with these estimators. 
Statistical studies for galaxy distribution have already shown that estimated values of these moments are more likely to have 
a lower value then its mean (Szapudi & Colombi 1996). This will mean that the probability distribution of these estimators 
are skewed. The skewness associated with these estimators can also be computed using the procedure outlined above and we 
hope to present a detailed analysis elsewhere; see Szapudi et al. (2000) for comments in a similar vein. 

The results presented above are valid for only one particular patch of the sky for computation of cumulants and a single 
pair of patches for measurements of cumulant correlators. However it is straightforward to generalize the above results for 
large number of patches or pairs of patches which will help to increase the signal to noise ratio (see Schneider et al. 1998). 

Detailed expressions for many-point cumulant correlators can also be obtained using the formalism developed here. The 
general expression we have presented in the text are valid for an arbitrary number of points and do not depend directly on 
the filter function used. In deriving the above results we have assumed that the two patches have same size and hence the 
same variance, but they might contain different number of galaxies, i.e. m and n 2 . It is not difficult to generalize the results 
to the case when the variance in these two patches are also different. We aim to present numerical results for specific filters 
which arise from above results and their comparison against simulated noisy convergence maps in the near future. 
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